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d 要 : 宁夏 沿 黄 城 市 带 是 宁夏 经 济 社会 发 展 的 精华 地 带 , 探 索 该 区 域 的 生态 环境 问题 对 宁夏 沿 黄 城 市 带 建 设 意 


义 重大 。 本 文 基于 2002 一 2018 年 的 MODIS 遥感 数据 ,利用 Mann-Kendall 非 参 数 检验 、Sen’s 斜率 估计 、 景 观 格局 分 


析 、 偏 相关 分 析 以 及 残 差 趋势 分 析 方 法 ,研究 了 该 区 域 植 被 覆盖 状况 的 时 空 演变 过 程 及 趋势 。 结 果 表 明 :2002 一 
2018 年 宁夏 沿 黄 城市 带 约 44.91% 区 域 的 植被 覆盖 度 显 著 增 加 , 约 1.31% 区 域 的 植被 覆盖 度 显 著 减 少 , 约 53. 7896 


区 域 的 植被 覆盖 度 变 化 不 显著 ;宁夏 治 黄 城市 带 植被 覆盖 度 变化 与 气温 具有 较 强 的 负 相关 性 ,与 降水 具有 较 强 的 


正 相 关 性 ;宁夏 沿 黄 城市 带 81. 4096 的 区 域 人 为 因素 对 植被 覆盖 度 多 年 变化 趋势 的 贡献 率 大 于 50% ,有 约 24% 的 
区 域 气候 因素 对 植被 覆盖 度 多 年 变化 趋势 的 贡献 率 大 于 50% 。 


关键 词 : 植被 覆盖 度 ; 时 空 演 变 ; 驱动 力 ; Mann-Kendall 非 参数 检验 ; Sen's 斜率 估计 ; 沿 黄 城 市 带 ; 宁夏 


宁夏 沿 黄 城市 带 是 以 地 缘 相 近 交通 便利 经济 
关联 度 较 高 的 银川 市 为 中 心 ,石嘴山 、 匡 忠 、 中 卫 3 
个 地 级 市 为 主干 形成 的 城镇 集合 体 ,是 宁夏 城镇 分 
布 最 密集 、 生 产 力 要 素 最 活路 的 区 域 ,也 是 宁夏 经 济 
社会 发 展 的 精华 地 带 。 宁 夏 沿 黄 城 市 带 以 宁夏 全 区 
43% 的 土地 面积 创造 了 该 区 90% 以 上 的 GDP 和 财 
政 收 入 ,承载 着 引领 宁夏 工业 化 和 城市 化 建设 的 重 
任 。 然 而 ,工业 化 和 城镇 化 的 迅速 发 展 ,给 社会 禹 来 
巨大 效益 的 同时 ,也 可 能 会 引发 诸多 生态 环境 问题 ， 
例如 大 气 污染 .水体 污 染 、 绿 地 覆盖 度 减 少 、 湖 泊 面 
积 缩小 等 。 定 量 评价 与 深入 分 析 沿 黄 城 市 带 工业 化 
和 城镇 化 已 经 或 引发 的 生态 环境 问题 ,对 于 宁夏 沿 
黄 城 市 带 建 设 具 有 重大 意义 。 植 被 作为 陆地 生态 系 
统 的 主体 ,是 连接 土壤 、 大 气 和 水 分 等 要 素 的 自然 纽 
带 , 因 此 成 为 陆地 生态 环境 的 重要 “指示 器 ”之 
一 0-31。 植 被 覆盖 度 (fractional vegetation. cover, 
FVC) 定 义 为 植被 垂直 投影 面积 与 地 表面 积 的 比值 ， 
通常 用 百分比 (% ) 来 表示 ,是 进行 区 域 生态 系统 评 
估 的 重要 定量 指标 "*““ 。 本 文 针对 宁夏 沿 黄 城市 带 
植被 覆盖 度 的 时 空 演变 过 程 和 趋势 开展 研究 ,期 望 
从 植被 覆盖 度 的 角度 ,透视 宁夏 沿 黄 城市 带 工业 化 
和 城镇 化 已 经 或 引发 的 生态 环境 问题 。 遥 感 技术 在 
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监测 植被 覆盖 度 方 面具 有 多 光谱 、 多 分 辨 率 、 多 时 
相 、 多 角度 .快速 成 像 等 优势 "… ,已 经 成 为 监测 全 
球 及 区 域 植 被 覆盖 变化 的 最 主要 手段 。 目 前 ,已 形 
成 多 种 植被 覆盖 度 遥 感 计算 方法 ,这 些 方法 大 致 可 
分 为 :Q 物理 模型 法 :如 一 种 利用 中 等 分 辨 率 与 高 
分 辩 率 卫星 图 像 相 结合 建立 植被 覆盖 度 提 取 模 型 ， 
以 获取 大 范围 植被 覆盖 度 的 方法 ” ;动态 贝 叶 斯 网 
络 方法 ,将 辐射 传输 模型 与 作物 生长 统计 模型 相 结 
合 , 从 中 分 辩 率 成 像 光 谱 仪 (MODIS ) 反 射 率 数据 中 
提取 植被 覆盖 度 ,综合 利用 遥感 数据 和 作物 生长 特 
VERO (SLE I ; 叶 面 积 指数 (leaf area index,LAI) 模 型 
反 演 法 ,可 准确 估算 出 乔 、 灌 、 草 覆盖 度 "-" © 
回归 模型 法 : 如 利用 归 一 化 植被 指数 ( normalized 
difference vegetation index, NDVI) 和 实测 的 FVC T5 
E REL rg 7 ;利用 实测 的 地 面 覆 盖 度 数据 与 
SAVI (soil adjusted vegetation index) 建立 了 锡 林 浩 
特地 区 植被 覆盖 度 的 线性 回归 模型 ;选取 了 包括 
DVI (simple difference vegetation indices)  NDVI、SA- 
VI SE 20 种 植被 指数 ,利用 回归 模型 估算 植被 覆盖 
E, O 混合 像 元 分 解法 :如 利用 像 元 二 分 模型 
估算 植被 覆盖 度 …“ ” ;使 用 亚 像 元 分 解 模型 进行 杆 
被 覆盖 度 的 动态 变化 分 析 "” ;采用 干枯 燃料 指数 


基金 项 目 : 宁夏 回族 自治 区 重点 研发 计划 一 般 项 目 (2018BEG03069 ) ;国家 自然 科学 基金 项 目 (41871338 ) 资助 


作者 简介 : 吴 加 敏 (1965 — ) , 男 , 地 图 制图 学 与 地 理 信 息 工 程 博 了 


E-mail; 188406255@ qq. com 
通讯 作者 : PEG. E-mail; 18336831802€ 163. com 


F ,教授 级 高 级 工程 师 ,研究 方向 为 资源 环境 遥感 应 用 技术 . 


http :// azr. xjegi. com 


Chinai AEHL 
ChinaXivâ fERBT 


吴 加 敏 等 :宁夏 沿 黄 城市 带 植被 覆盖 时 空 演 变 及 其 驱动 力 分 析 


( dead fule index, DFI) 构建 NDVI-DFI 像 元 三 分 模 
型 ,估算 了 锡林郭勒 草原 的 光合 植被 覆盖 度 和 非 光 
合 植被 覆盖 度 '" ;采用 面向 对 象 多 端 元 混 解 模型 ， 
有 效 地 减少 了 计算 量 和 混合 像 元 的 端 元 变化 ,提高 
了 植被 覆盖 度 反 演 精度 。@ 机 器 学 习 法 :如 利 
用 神经 网 络 和 PROSAIL 辐射 传输 模型 来 反 演 植被 
和 覆盖 度 , 并 生成 了 CYCLOPES 3H zc RE j^ 4 UU ; DH 
支持 向 量 机 (support vector machine ,SVM ) 方法 来 识 
别 半 干旱 区 的 植被 类 型 ,在 SVM 结果 上 进一步 估算 
研究 区 域 的 FVC ;基于 IKONOS 遥感 影像 ,利用 
决策 树 分 类 方法 估算 乔木 植被 覆盖 度 ,分 类 精度 达 
97.395 91. 

本 文 的 主要 目的 在 于 揭示 宁夏 沿 黄 城市 带 
FVC 的 时 空 演变 过 程 和 趋势 ,考虑 到 长 时 间 序 列 分 
析 对 遥感 数据 的 一 致 性 要 求 .研究 区 的 地 表 和 覆盖 类 
型 特点 (主要 是 草地 和 耕地 ) ,以 及 技术 方法 的 鲁 棒 
性 ,可 行 性 等 ,基于 LAI 反 演 FVC 的 物理 模型 。 


1 数据 资料 


1.1 研究 区 概况 

图 1 为 宁夏 沿 黄 城市 带 的 空间 分 布 ,图 中 的 土 
地 和 履 盖 类 型 为 2017 年 的 MODIS (or moderate resolu- 
tion imaging spectroradiometer ) 土地 覆盖 类 型 产品 
MCDI2QI ,分 类 模式 为 IGBP 全 球 植被 分 类 ,该 土地 
覆盖 类 型 的 空间 分 辨 率 是 500 m。 宁 夏 沿 黄 城市 带 
国土 面积 约 为 4.16 x 10* km? ,地 理 坐 标 104?22' ~ 
107°50'E、36°43' ~39°22'N, 属 温带 大 陆 性 半 湿 润 半 
干旱 气候 ,年 均 气 温 5 ~9 "C ,年 降水 180 ~680 mm, 
该 区 域 的 地 表 覆 盖 类 型 主要 是 草地 和 耕地 , 占 整个 
研究 区 的 9396 以 上 ,另外 还 有 少量 的 灌 从 ,城市 和 
建成 区 以 及 稀 玻 植被 。 
1.2 遥感 数据 

本 研究 主要 采用 MODIS 叶 面 积 指数 (LAI) 产 
品 ` 地 表 覆 盖 类 型 产品 ,太阳 天 项 角 数 据 产 品 以 及 地 
表 反 射 率 产 品 。 其 中 LAI 产品 具体 为 MCD15A2 ,时 
间 分 辨 率 为 8 d, 空 间 分 辩 率 为 500 m, ,时 间 范 围 为 
2002—2018 年 ,影像 合计 759 景 ,空间 轨道 范围 为 
h26v05 。 地 表 覆 盖 类 型 产品 具体 为 MCD12Q1 ,时 间 
分 辨 率 为 1 a, 空 间 分 辩 率 为 500 mm ,时间 为 2017 
年 ,影像 合计 1 景 。 太 阳 天 顶 角 数据 产品 具体 为 
MODO9A1 ,时 间 分 辨 率 为 8 d, 空 间 分 辨 率 为 300 m, 
时 间 范 围 为 2002—2018 年 ,影像 合计 759 景 。 地 表 
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图 1 宁夏 沿 黄 城市 带 空 间 分 布 以 及 主要 地 表 和 覆盖 类 型 


Fig.1 Ningxia along the yellow city belt distribution 


反射 率 产品 为 MOD09Al ,时间 分 别 为 MODOO AT. 
A2018177, MODO9Al1.  A2018193, MODO9Al. 
42018209, MODO9A1. A2018225 共 4 景 。 首 先 利用 
MODIS Reprojection Tools ( MRT) 工具 对 LAI 数据 、 
太阳 天 项 角 数 据 .地 表 履 善 类 型 产品 ;以 及 地 表 反 射 
率 数 据 进行 坐标 转换 、 格 式 转换 .投影 转换 等 一 系列 
预 处 理 。 通 过 ENVI/IDL, 利用 研究 区 的 矢量 数据 
进行 裁剪 得 到 研究 区 LAI 数据 地表 覆盖 类 型 产品 、 
太阳 天 项 角 以 及 地 表 反 射 率 数据 。 


2 研究 方法 


2.1 植被 覆盖 度 计算 
基于 辐射 传输 模型 理论 ,植被 覆盖 度 与 植被 叶 
面积 指数 密切 相关 ,存在 一 个 经 典 物理 意义 上 的 转 
BAR”, 
FVC =1 — el -6(0) x LAVcost] (1) 
RP: G(9) 为 叶片 投影 因子 , 取 值 为 0. 355 ; 9 为 
太阳 天 项 角 。 
FVC zi ze DAV (26088) (2) 
式 中 :FVC 的 取 值 范围 为 0 ~1, 通 常情 况 下 FVC < 
20% 为 低 植被 覆盖 区 ;20% < FVC <40% 为 中 低 植 
被 覆盖 区 ;40% < FVC < 60% 为 中 植被 覆盖 区 ; 


| 


于 oH 


6096 « FVC <80% 为 中 高 植被 覆盖 区 ; FVC > 8096 
为 高 植被 覆盖 区 。 

此 外 ,为 减少 大 气 的 云 、. 颗 粒 阴影、 视角 以 及 太 
阳 高 度 角 等 对 LAT 数据 的 影响 ,本 文 根 据 具体 的 
人 研究 尺度 ,采用 最 大 值 合 成 法 (maximum value com- 
posit, MVC) ,对 8 d 尺度 的 LAI 进一步 处 理 ,其 表达 
式 如 下 : 


d 
LAT = MaxLAI, (3) 


式 中 :LAT 是 第 i 年 (月 )LAI 最 大 合成 值 , 即 LAT 的 
最 高 值 ;LAL 为 各 像 元 8 d 合成 最 大 LAI 值 ;d 为 第 
年 (月 )8 d EI LAIL 影像 总 数 。 

2.2 时 间 序 列 分 析 方 法 

2.2.1  Mann-Kendall 非 参 数 检 验 Mann-Kendall 
非 参数 检验 方法 ”是 目前 被 广泛 应 用 的 趋势 分 析 
方法 ,能 很 好 地 揭示 时 间 序 列 的 趋势 变化 ,其 优点 在 
于 不 需要 遵从 一 定 的 正 态 分 布 ,也 不 受 少数 异常 值 
的 干扰 。 设 某 时 间 序 列 为 1X,| (i=1,2,…,n), 时 
EEJ X} 的 对 侦 数 : 


n-l n 
S= 22 Sen CX; -X,) 


+1,(X,-X;)>0 
sgn( X, - X;) - 0,(X,—X,) =0 

-1,(X, -X,) «0 

式 中 :5 为 正 态 分 布 ,均值 为 0 ,方差 Var(S) 2n(n- 

1) (2n * 5)/18, ?5 n > 10 时 , Mann-Kendall 统计 


Æ U: 


(4) 


La S>0 
/ Var( S) i 
U= 0, S=0 (5) 
S+l 
—————, S«0 
/ Var( S) 


M-K 统计 量 0 的 取 值 范围 为 ( -% , +o), 4H 
果 1071> Do， 即 “置信 水 平 上 ,时 间 序 列 数据 存在 
明显 的 上 升 或 下 降 趋 势 。U > 0 时 ,表示 时 间 序 列 
(XI AER; U «0 时 ,表示 时 间 序 列 |X,} 为 下 
降 趋势 。 

根据 显著 性 水 平 以 及 统计 量 0 将 变化 趋势 分 
为 5 个 等 级 ,分 级 结果 见 表 1。 
2.2.2 Sen's 斜率 估计 Mann-Kendall 算法 可 以 实 
现时 间 序 列 的 趋势 显著 性 检验 ,但 无 法 得 到 时 间 序 
列 的 变化 率 。Sen’'s 斜 率 估计 对 Mann-Kendall 算 法 


"umm 
Æ HO- 
21 


R1 显著 性 检验 结果 分 级 


Tab.1 Classification of significance test results 


变化 趋势 统计 量 U 显著 性 水 平 
极 显 著 减 少 U «0 a «0.01 
一 般 显著 减少 0.01 <a «0.05 
基本 不 变 o >0.05 
一 般 显 著 增 加 U»0 0.01 <a «0.05 
极 显著 增加 o «0.01 


进行 了 补充 ,能 够 准确 计算 时 间 序列 的 变化 率 B 
当 B>0 时 ,序列 旦 上升 趋 势 ; 当 B 20 时 ,表示 序列 
变化 趋势 不 明显 ; 当 有 <0 时 ,序列 时 下 降 趋势 。 

对 于 时 间 序 列 x; = (xi x, 7x, ) ,Sen's SER B 
的 计算 公式 为 


X; 一 多 


B = Median ( visi (6) 


J=! 
式 中 ;Median Jy P fft 8c, BD AEA E 0: 3 c E 
从 小 到 大 的 顺序 排列 ,如 果 序 列 对 是 奇数 , 则 取 最 中 
间 的 一 个 数 作为 Sen's 斜率 ,如 果 序 列 对 是 偶数 , 则 
取 最 中 间 的 两 个 数 的 平均 数 作为 Sen's 斜率 ;Y 表 
示 任 意 。 

2.3. ，” 偏 相关 分 析 以 及 残 差 趋势 分 析 

2.3.1. 偏 相关 分 析 ” 偏 相关 分 析 是 指 两 个 变量 同 
时 与 第 3 个 变量 相关 时 ,将 第 3 个 变量 的 影响 剔除 ， 
只 分 析 另 外 两 个 变量 间 相关 程度 的 过 程 * ,该 方法 
在 一 定 程度 上 能 够 减弱 因 某 一 因素 的 变化 导致 其 他 
因素 与 因 变量 的 关系 的 影响 。 本 文 基于 像 元 计算 
FVC 与 气候 指标 的 偏 相 关系 数 ,首先 计算 FVC 与 年 
平均 气温 和 年 降水 的 简单 相关 系数 ,进而 得 到 偏 相 
关系 数 。 计 算 公式 如 下 : 

(0 RQ-R.OR, 
VO-RJO-R) 


式 中 ;R,, .为 固定 变量 后 变量 x 与 y 的 偏 相关 系数 ; 
Ro Ra RIIN x 与 y,x 与 z,y 与 z 之 间 的 相关 系 
2.3.2. 残 差 趋势 分 析 气候 因素 以 及 人 为 因素 是 
影响 FVC 变化 的 两 大 因素 。 其 中 气温 、 降 水 是 气候 
因素 中 的 两 大 主要 因素 。 鉴 于 其 他 气候 因素 对 
FVC 变化 的 影响 相对 较 小 ,笔者 在 气象 因素 中 仅 考 
虑 气温 .降水 两 大 因素 。 利 用 残 差 趋势 分 析 法 将 每 
个 像 元 的 FVC 与 气候 因子 (年 平均 气温 .年 降水 ) 做 
回归 分 析 ,进而 得 到 每 个 像 元 的 FVC 预测 值 , 该 预 
测 值 为 单个 气候 因子 对 FVC 的 影响 ,以 偏 相 关系 数 


(7) 


x yz 
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的 大 小 作为 权 值 进行 加 权 平 均 , 获 得 最 终 FVC 的 预 
测 值 , 即 气候 因素 对 FVC 的 贡献 值 。 将 每 个 像 元 的 
FVC 真 值 减 去 FVC 预测 值 可 到 人 类 活动 对 FVC 的 
影响 值 , 即 人 为 因素 对 FVC 的 贡献 值 。 具 体 为 : 


FVC uas = W pred temp X FVC pred temp + W predi pree x 
FVC pred prec ( 8 ) 
式 中 下 = R pyc temp. prec wW 
T : pred temp ^ R 十 R , pred. prec 
FVC temp. prec FVC prec. temp 
— Rpve prec. temp FVC FVC VAN 
= ; pred temp » pred. prec JJ 


Rpvye temp. pree + BEVC. prec_temp 
别 为 FVC 与 年 平均 气温 .年 降水 线性 回归 预测 值 。 

FV Chaman y = FVC ea -FVC imaa (9) 

IRP : FV Cpm ICEP t FAHARI FVC 的 

贡献 值 ; FVC。 为 像 元 在 第 上 年 的 真实 值 ; 

FVC aima ,为 像 元 在 第 上 年 气候 因素 对 FVC 的 贡献 


true £ 


值 。 
dE VG... dFVC,, a. dE VEL 
= + 
dt dt dt 
= C( climate) + C( human) (10) 
m 中 : dFVC ue” dt ` dFVC LZ dt ` dFVC puman” dt 分 别 
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图 2  LAI scii FVC-NDVI Ý FVC 交叉 验证 
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表示 FVC 多 年 变化 趋势 .气候 因素 、 人 为 因素 对 
FVC 多 年 变化 趋势 的 贡献 。 

进一步 可 以 求 得 气候 以 及 人 为 因素 对 FVC 多 
年 变化 趋势 的 贡献 率 p(Y) 为 : 


_ C) 
PX) - wc. 7di 


3 结果 与 讨论 


3.1 植被 覆盖 度 计算 结果 交叉 验证 

利用 像 元 二 分 模型 计算 FVC ,根据 沿 黄 城市 带 
植被 覆盖 的 实际 情况 以 及 植被 指数 NDVI 灰 度 分 布 
频率 累计 表 , 选 择 99.5% 和 0. 5% 的 置信 区 间 内 上 
下 限 累 计 频 率 值 作为 纯 植 被 覆盖 时 的 NDVI 和 纯 土 
壤 地 表 的 NDVI。 从 像 元 二 分 模型 和 LAI 模型 中 等 
间隔 提取 200 个 像 元 对 ,相关 性 分 析 如 图 2 所 示 。 
结果 表明 两 者 的 相关 系数 分 别 为 :0. 848 (177), 
0. 847 (193 ) ,0. 766 (209 ) ,0.775(225 ) , 均 在 显著 性 
水 平 P <0.01 的 情况 下 显著 相关 ;其 次 斜率 分 别 为 : 
1.046(177) ,0.974(193) ,0.923(209) ,0.971(225) , 


x 10096 (11) 
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Fig.2 LAI retrieve FVC-NDVI retrieve FVC cross validation 
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到 3 2002—2018 年 宁夏 沿 黄 城市 带 植被 覆盖 度 分 布 


Fig.3 Vegetation coverage distribution along the yellow city belt in Ningxia from 2002 to 2018 


均 接 近 于 y =x。 因 此 ,本 文 基于 LAIL 计算 的 FVC H. 
有 较 大 的 可 靠 性 。 
3.2 植被 覆盖 度 的 空间 分 布 格局 

图 3 为 2002—2018 年 (间隔 4 a 显示) 宁夏 沿 
黄 城市 带 植被 覆盖 度 空间 分 布 。 可 以 看 出 ,石嘴山 
市 南部 植被 覆盖 度 较 高 ,北部 偏 低 ;银川 市 植被 覆盖 
度 最 高 , 且 主 要 集中 在 北部 ;吴忠 市 植被 覆盖 度 整 体 
偏 低 , 较 高 区 域 主 要 集中 在 西北 部 ;中 卫 市 植被 覆盖 
度 整体 偏 低 , 较 高 区 域 主要 集中 在 西北 部 。 其 次 ， 
2002—2006 年 植被 覆盖 度 最 低 。 

图 4 为 宁夏 沿 黄 城市 带 2002—2018 年 (间隔 
4 a) 各 等 级 植被 覆盖 面积 占 比 。 从 图 中 可 以 看 出 ， 
5 个 年 份 植被 覆盖 度 低 的 区 域 面积 占 比较 高 。 其 中 
低 植被 覆盖 区 面积 占 比 在 2006 年 出 现 了 极 显 著 增 
加 ,达到 了 44. 6196 ;中 低 植被 覆盖 区 面积 占 比 出 现 极 
显著 减少 ,减少 至 34. 66% ;其 余 4 个 年 份 面积 占 比 相 
对 稳定 ;中 .中 高 ,高 植被 覆盖 区 面积 占 比 在 这 5 a 间 
呈 增 加 的 趋势 ,分 别 由 2002 年 的 14. 94% ,8. 0496 、 
2.79% 增 加 至 2018 年 的 31.18% ,14. 9296 .7.84% 。 
3.3 植被 覆盖 度 的 时 间 变 化 过 程 

MODIS LAI 产品 起 始 时 间 为 2002 年 7 月 4 日 ， 


低 植 被 覆盖 。 多 绝 中 低 植 被 覆盖 ES 中 植被 覆盖 
中 高 植被 覆盖 高 植被 覆盖 
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图 4 2002—2018 年 各 等 级 植被 覆盖 度 面积 占 比 
Fig.4  Proportion of vegetation coverage of each 
grade from 2002 to 2018 


故 2002 4E 6 月 无 相应 的 数据 。 由 图 5 可 见 ,植被 覆 
HEE 5 一 9 月 呈现 出 先 上 升 后 下 降 的 变化 趋势 ， 
5 一 8 月 植被 覆盖 度 有 所 上 升 , 且 8 月 植被 覆盖 度 达 
到 最 大 ,9 月 植被 覆盖 度 显著 下 降 。 其 次 ,5 个 年 份 
中 ,2018 年 各 月 份 植被 覆盖 度 均 最 大 ; 而 2006 年 各 
月 份 植被 覆盖 度 均 最 小 。 

由 图 6 可 知 ,2002 一 2018 年 年 均 植 被 覆盖 度 总 
体 呈 线性 增长 态势 (0.78% « a^ ,P «0.01) ,最 低 
值 27.68% 发 后 在 2005 年 ,最 大 值 46. 1296 发 生 在 
20184£ , HP, 2002—20054£ © F EE d 38, H 
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图 5 2002—2018 年 生长 季 ( 间 隔 4 a) FVC 年 内 变化 
Fig.5 Annual changes of FVC during the growth season 


from 2002 to 2018 (4 years interval) 
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图 6 2002—2018 年 年 均值 FVC 变化 
Fig.6 Annual mean FVC change from 2002 to 2018 


33.81% 下降 至 27. 68% ;2005 一 2007 一 2008 年 呈现 
先 增 8.52% 后 降 4.52% ;2008 一 2012 一 2015 年 呈 先 
增 10. 8396 后 降 7. 24% ; 2015 一 2018 年 呈 上 升 趋 
势 , 由 35.27% 上升 至 46.12% 。 
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3.4 植被 覆盖 度 的 演变 趋势 

图 7 为 2002—2018 年 各 等 级 植被 覆盖 度 时 间 
序列 变化 。 总 体 看 来 ,2002 一 2018 年 低 植被 覆盖 像 
元 数 呈 减少 趋势 ( -2 748 . a ,P<0.05); 中 低 植 
被 覆盖 像 元 数 增幅 极 小 且 极 不 显著 ;中 植被 覆盖 区 、 
中 高 植被 覆盖 区 ,高 植被 覆盖 区 FVC 像 元 数 均 呈 增 
加 趋势 。 其 中 ,中 高 植被 覆盖 区 高 植被 覆盖 区 像 元 
数 增加 趋势 极 显著 (P <0.01)。 

图 8 是 利用 Mann-Kendall 非 参 数 检 验 和 Sen's 
斜率 分 析 得 到 的 宁夏 沿 黄 城市 带 植被 覆盖 度 变 化 趋 
势 及 变化 率 , 表 2 统计 了 相应 的 面积 变化 及 其 比例 。 
研究 表明 ,2002 一 2018 年 宁夏 沿 黄 城市 带 植 被 覆盖 
度 显 著 增加 的 区 域 占 总 面积 的 44. 9196 , 主要 分 布 
在 清水 河和 黄河 沿岸 ,以 及 石嘴山 和 中 卫 市 的 东南 
部 区 域 ; 植 被 覆盖 度 显著 减少 的 区 域 占 总 面积 的 
1.31% ,主要 分 布 在 银川 .吴忠 和 中 卫 市 建成 区 及 其 
附近 ,主要 是 由 于 城市 建筑 用 地 扩张 引起 的 ;植被 覆 
盖 度 变化 不 显著 的 区 域 占 总 面积 的 53. 78960 , Hc 
中 ,植被 覆盖 度 一 般 显著 增加 和 极 显著 增加 的 面积 
分 别 占 研 究 区 总 面积 的 19. 37% 和 25. 54% 。 相 比 
之 下 ,植被 覆盖 度 减 少 的 地 区 面积 占 比 相对 较 小 。 
其 中 ,植被 覆盖 度 一 般 显著 减少 和 极 显 著 减 少 的 面 
积分 别 占 研究 区 总 面积 的 0.65% 和 0. 66% 。FVC 
变化 率 小 于 -0.000 4 的 区 域 占 总 面积 的 6.78% ; 
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图 7 2002—2018 年 不 同 植被 覆盖 像 元 数 变 化 


Fig.7 Variation of pixel number of different vegetation cover from 2002 to 2018 
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8 2002—2018 宁夏 沿 黄 城 市 带 FVC 变化 趋势 及 变化 率 
Fig.8 Change trend and change rate of FVC along yellow city belt in Ningxia from 2002 to 2018 


R2 2002—2018 年 宁夏 沿 黄 城市 带 Sen's +M - K 面积 占 比 变化 


Tab.2 Area change of Sen's +M -K in the yellow belt of ningxia from 2012 to 2018 / 96 
( 0.048, -0.01] ( -0.01, -0.0004 ] ( —0.0004 ,0. 0004 ] (0. 0004 ,0. 01 ] (0.01,0. 062] 合计 
极 显著 减少 0.56 0.10 0 0 0 0. 66 
一 般 显著 减少 0.31 0.34 0 0 0 0. 65 
基本 不 变 0.20 5.27 2.79 43.39 2.13 53.78 
一 般 显著 增加 0 0 0 12. 12 7.25 19.37 
极 显著 增加 0 0 0 5.11 20.43 25. 54 
合计 1.07 5.71 2.79 60. 62 29.81 100. 00 


变化 率 介 于 -0.000 4 ~0.000 4 之 间 的 区 域 占 研究 
区 总 面积 的 2.79% ;变化 率 介 于 0.000 4~0.01 的 
区 域 面积 占 研究 区 总 面积 的 60. 62% ;变化 率 介 于 
0.01 ~0. 062 的 区 域 占 研究 区 总 面积 的 29. 81% 。 
总 之 ,从 变化 趋势 以 及 变化 率 的 角度 均 可 看 出 ,宁夏 
治 黄 城市 带 整 体 FVC 增加 趋势 较为 明显 。 


4 讨论 


4.1 植被 覆盖 度 变 化 与 气候 因素 的 相关 性 

图 9 为 宁夏 沿 黄 城 市 带 年 最 大 FVC 与 年 平均 
气温 和 年 降水 的 偏 相关 系数 空间 分 布 。 由 年 最 大 
FVC 与 年 平均 气温 的 空间 相关 性 分 析 可 知 , 正 相关 
的 区 域 主要 分 布 在 黄河 以 及 清水 河沿 岸 周边 ,土地 
利用 类 型 以 耕地 为 主 ;呈现 负 相关 和 较 弱 相关 的 区 
域 主要 集中 在 城市 和 建成 区 周边 ,以 及 中 卫 市 中 部 
和 东南 部 区 域 ,土地 利用 类 型 以 草原 为 主 。 由 年 最 
大 FVC 与 年 降水 空间 相关 性 分 析 可 知 ,与 气温 的 空 
间 相关 性 相反 , 正 相关 区 域 主要 集中 在 中 卫 市 \ 匡 上 忠 


市 以 及 银川 市 的 西北 部 和 东南 部 ,主要 土地 利用 类 
型 为 草地 , 少 部 分 为 耕地 。 呈 现 负 相关 以 及 相关 性 
较 弱 的 区 域 主要 集中 在 城市 和 建成 区 附近 。 整 体 而 
言 , 与 年 均 气 温 相 比 ,年 最 大 FVC 与 年 降水 的 相关 
性 更 强 。 
4.2 ”气候 及 人 为 因素 对 植被 覆盖 度 趋势 的 贡献 
植被 覆盖 度 变 化 是 气候 因素 和 人 为 因素 共同 作 
用 的 结果 ,分 析 气 候 因素 和 人 为 因素 对 宁夏 沿 黄 城 
市 带 植 被 覆盖 度 的 变化 趋势 贡献 尤为 重要 。 图 10 
为 气候 及 人 为 因素 对 多 年 FVC 变化 趋势 贡献 率 的 
空间 分 布 。 从 图 10a 可 知 ,贡献 率 小 于 11% 的 区 域 
主要 集中 在 城市 建成 区 及 其 周围 , 面积 占 比 约 
6.8% ,其 中 石嘴山 市 较为 集中 ;贡献 率 介 于 11% ~ 
31% 之 间 的 区 域 主要 集中 在 黄河 沿岸 以 及 清水 河沿 
岸 ,面积 占 比 约 37.7% 。 贡 献 率 大 于 50% 的 区 域 呈 
零散 分 布 ,面积 占 比 约 24% ,在 中 卫 市 .吴忠 市 较为 
集中 。 从 图 10b 可 知 ,6. 2% 的 区 域 贡 献 率 小 于 
3196 ; 12. 496 的 区 域 贡 献 率 介 于 31% ~50% 之 间 ; 
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图 9 宁夏 沿 黄 城市 带 年 最 大 值 FVC 与 年 平均 气温 和 年 降水 的 偏 相 关系 数 空间 分 布 


Spatial distribution of partial correlation coefficient between annual maximum FVC and annual mean temperature 


and annual precipitation in the yellow city belt of Ningxia 
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图 10 2002—2018 年 宁夏 沿 黄 城市 带 气 候 因 素 和 人 为 因素 对 FVC 变化 趋势 的 贡献 率 
Fig. 10 Contribution rate of climate factors and human factors to FVC variation trend in the yellow city belt of Ningxia 


from 2002 to 2018 


81.496 的 区 域 贡 献 率 大 于 50%。 相 比 于 气候 因素 ， 
人 为 因素 对 FVC 的 贡献 率 整 体 较 高 。 


5 结论 


(1) Mann-Kendall 非 参 数 检验 表明 ,2002 一 
2018 年 宁夏 沿 黄 城市 带 约 44.91% 区 域 的 植被 覆盖 
度 显著 增加 , 约 1.31% 区 域 的 植被 覆盖 度 显 著 减 
少 , 约 53.78% 区 域 的 植被 覆盖 度 变 化 不 显著 。 


(2) Sen's 斜率 分 析 表 明 ,2002 一 2018 年 宁夏 
沿 黄 城 市 带 约 60. 61% 区 域 的 FVC 变化 率 介 于 
0.000 4 ~0.01 之 间 , 约 29.81% 区 域 的 FVC 变化 率 
介 于 0.01 ~0.062 之 间 , 约 2.80% 区 域 的 FVC 变化 
率 介 于 -0.000 4 ~0.000 4 之 间 , 约 6. 7896 区域 的 
FVC 变化 率 小 于 -0.000 4, 

(3) 宁夏 沿 黄 城 市 带 FVC 与 年 平均 气温 叶 显 
著 负 相 关 ,与 年 降水 呈 显 著 正 相关 ,说 明 气候 因子 是 
影响 宁夏 沿 黄 城市 带 FVC 的 重要 因素 之 一 。 


v:202006.00190v1 


E mM 


chinaX 


TOS 区 


(4) 植被 覆盖 度 变 化 是 气候 因素 和 人 为 因素 共 


同 作 用 的 结果 ,81. 496 的 区 域 人 为 因素 对 植被 覆盖 
度 多 年 变化 趋势 的 贡献 率 大 于 50% ,有 约 24% 的 区 


域 所 


候 因素 对 植被 覆盖 度 多 年 变化 趋势 的 贡献 率 大 


于 50% ,人 为 因素 已 成 为 宁夏 沿 黄 城市 带 杆 被覆 疼 
度 变 化 的 主导 因素 。 


总 之 ,2002 一 2018 年 宁夏 沿 黄 城 市 带 植 被 覆盖 


度 总 体 呈 增 大 趋势 ,其 中 植被 覆盖 度 显著 增加 的 区 
域 主要 分 布 在 清水 河和 黄河 沿岸 以 及 石嘴山 和 中 卫 


市 的 东南 部 ;显著 减少 区 域 主要 分 布 在 银川 .吴忠 和 
中 卫 市 建成 区 及 其 附近 ,主要 是 由 于 城市 建筑 用 地 
扩张 引起 的 。 
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Spatiotemporal evolution and driving force analysis of fractional vegetation 
coverage over the urban belt along the Yellow River in Ningxia 


WU Jia-min', CAI Chuang-chuang , SUN Hao’, YAO Jian-hua! , 
CHEN Wei’, GU Ji-sheng', JIANG Jin-bao’ 
(1. The Remote Sensing Surveying and Mapping Institute of the Ningxia Hui Autonomous Region/ Remote Sensing Center of Ningxia , 
Yinchuan 750021 , Ningxia , China ; 
2. College of Geoscience and Surveying Engineering , China. University of Mining and Technology ( Beijing) , Beijing 100083 , China ) 


Abstract: The urban belt along the Yellow River in Ningxia is the essence of Ningxia's economic and social de- 
velopment. It is significant to study the environmental issues in these areas. In this study , we investigated the spatio- 
temporal evolution and trend of fractional vegetation coverage ( FVC) based on MODIS remote sensing data from 
2002 to 2018 using the Mann-Kendall non-parametric test, Sen's slope estimation, and landscape pattern analysis 
methods. The results showed that from 2002 to 2018 , FVC in 44. 9196 of the study area increased significantly , FVC 
in 1. 3196 of the area decreased ,and FVC did not change in 53. 7896 of the area. There was a strong negative corre- 
lation between vegetation coverage change and temperature and a strong positive correlation between vegetation cov- 
erage change and precipitation. The contribution of 81. 4096 of regional man-made factors along the Yellow City Belt 
in Ningxia to the multi-year trend of FVC was greater than 5096 , and about 24% of the regional climate factors con- 
tributed more than 50% to the multi-year trend of FVC. 

Key words: fractional vegetation cover; spatiotemporal evolution; driving force; Mann-Kendall non-parametric 


test; Sen's slope estimator; urban belt along the Yellow River; Ningxia 


